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A  two-dimensional  model  is  developed  to  study  the  performance  of  a  planar  solid  oxide  fuel  cell  (SOFC) 
running  on  steam/methane  mixture.  The  model  considers  the  heat/mass  transfer,  electrochemical  reac¬ 
tions,  direct  internal  reforming  of  methane  (CH4),  and  water  gas  shift  reaction  in  an  SOFC.  It  is  found  that 
at  an  operating  potential  of  0.8  V,  the  upstream  and  downstream  of  SOFC  work  in  electrolysis  and  fuel 
cell  modes,  respectively.  At  the  open-circuit  voltage,  the  electricity  generated  by  the  downstream  part 
of  SOFC  is  completely  consumed  by  the  upstream  through  electrolysis,  which  is  contrary  to  our  common 
understanding  that  electrochemical  reactions  cease  under  the  open-circuit  conditions.  In  order  to  inhibit 
the  electrolytic  effect,  the  SOFC  can  be  operated  at  a  lower  potential  or  use  partially  pre-reformed  CH4 
as  the  fuel.  Increasing  the  inlet  gas  velocity  from  0.5 ms-1  to  5.0ms_1  does  not  reduce  the  electrolytic 
effect  but  decreases  the  SOFC  performance. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  can  directly  convert  the  chemical 
energy  of  a  fuel  into  electricity  through  electrochemical  reaction 
in  an  efficient  and  environmentally  friendly  way.  An  SOFC  usu¬ 
ally  employs  oxygen-ion  conducting  ceramics  as  electrolyte,  like 
yttrium-stabilized-zirconia  (YSZ)  [1-3].  As  the  electrolyte  materi¬ 
als  exhibit  moderate  oxygen  ion  conductivity  at  a  high  temperature, 
the  operating  temperature  of  SOFC  is  usually  high  (i.e.  873-1273  K). 
The  high  working  temperature  makes  SOFC  very  suitable  for  com¬ 
bined  heat  and  power  co-generation,  as  the  waste  heat  from  SOFC 
is  of  high  quality  and  can  be  recovered  by  integrating  the  SOFC 
stack  with  absorption  heat  pumps  or  other  thermodynamic  cycles 
to  achieve  a  higher  overall  efficiency  [4-6].  Thermodynamic  anal¬ 
yses  have  been  reported  that  the  overall  efficiency  of  SOFC/gas 
turbine  system  can  be  higher  than  70%  [7-9].  In  addition,  carbon 
monoxide  (CO),  a  poisonous  gas  for  low  temperature  fuel  cells  (i.e. 
proton  exchange  membrane  fuel  cell:  PEMFC)  [10],  does  not  poison 
the  anode  catalyst  of  SOFC.  Instead,  CO  can  be  used  as  a  fuel  in  an 
SOFC.  Thus  alternative  fuels,  like  methane,  methanol,  and  ethanol, 
can  be  used  in  SOFCs  for  power  generation  [11-21].  The  fuel  flexi¬ 
bility  makes  SOFCs  more  advantageous  than  hydrogen-fueled  fuel 
cells,  as  it  is  still  very  difficult  to  produce  and  store  hydrogen  effec¬ 
tively  and  economically.  Due  to  their  great  prospect  for  clean  power 
generation,  SOFCs  have  received  more  and  more  attention  and 
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extensive  research  works  have  been  conducted  in  recent  years  to 
improve  the  SOFC  efficiency  and  its  long  term  stability  [22-27]. 

Among  a  lot  of  alternative  fuels,  methane  (CH4)  is  extensively 
studied  for  SOFC  as  it  is  the  major  component  in  natural  gas  (about 
90%)  and  one  of  the  major  components  in  biogas  [28].  As  direct  elec¬ 
trochemical  oxidation  of  methane  is  still  difficult  in  SOFC,  methane 
is  usually  reformed  either  externally  or  internally.  External  reform¬ 
ing  requires  additional  fuel  processing  units  for  methane  steam 
reforming  and  water  gas  shift  (WGS)  reactions,  thus  increasing 
the  system’s  overall  cost  and  complexity  [29-35].  For  comparison, 
direct  internal  reforming  (DIR)  eliminates  the  need  of  an  external 
reformer  as  the  high  working  temperature  enables  DIR  reaction  as 
well  as  WGS  reaction  for  H2  production  [36-52].  Thus  the  SOFC  sys¬ 
tem  can  be  simpler  and  compact.  In  addition,  part  of  the  thermal 
energy  requirement  for  DIR  reaction  can  be  provided  by  heat  gen¬ 
eration  in  SOFCs,  thus  a  higher  overall  efficiency  can  be  obtained. 
However,  the  endothermic  DIR  reaction  and  exothermic  WGS  reac¬ 
tion  complicate  the  SOFC  temperature  field  and  more  importantly, 
both  DIR  and  WGS  reactions  require  steam  in  the  anode,  thus  dilut¬ 
ing  the  fuel  concentration  and  reducing  the  SOFC  performance. 
In  addition,  carbon  deposition  may  occur  in  a  hydrocarbon-fueled 
SOFC  with  DIR  and  WGS,  which  in  turn  can  deteriorate  the  SOFC 
performance  considerably  [53].  In  order  to  solve  the  problem  of  car¬ 
bon  deposition  and  long  term  stability,  extensive  research  works 
have  been  conducted  both  experimentally  and  theoretically.  In 
addition  to  the  development  of  new  anode  materials,  one  com¬ 
mon  way  of  tackling  the  carbon  deposition  problem  is  to  supply 
sufficient  amount  of  steam  in  the  anode.  Based  on  thermodynamic 
analysis  in  the  literature,  the  steam-carbon  ratio  of  no  less  than 


0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 0.09.069 


2028 


M.  Ni  /  Journal  of  Power  Sources  196(2011)  2027-2036 


Nomenclature 

Bg  Permeability  of  the  porous  electrodes  (m2) 

cp  heat  capacity  (J  kg-1  K-1 ) 

da  thickness  of  anode  ([Jim) 

dc  thickness  of  cathode  ( pm) 

De^m  effective  diffusion  coefficient  of  species  i  in  gas  mix¬ 

ture  (cm2  s-1 ) 

Di  k  Knudsen  diffusion  coefficient  of  i  (cm2  s-1 ) 

Df  j  binary  diffusion  coefficient  of  i  and  j  (cm2  s-1 ) 

E  equilibrium  potential  (V) 

E0  reversible  potential  at  standard  condition  (V) 

F  Faraday  constant  (9.6485  x  1 04  C  mol-1 ) 

HDir  heat  demand  for  direct  internal  reforming  of 

methane  (J  mol-1 ) 

HWgs  heat  generation  from  water  gas  shift  reaction 
(J  mol-1 ) 

J  current  density  (A  m-2 ) 

k  thermal  conductivity  (W  m-1  K-1 ) 

Mf  molecular  weight  of  species  i  (kg  mol-1 ) 

P  operating  pressure  (bar) 

Pj  partial  pressure  (bar)  of  species  i  at 

electrode-electrolyte  interface 
Ract  resistivity  due  to  electrochemical  reaction  ( Q  m2 ) 

Rdir  reaction  rate  of  direct  internal  reforming  of  methane 

(molm-3  s-1) 

Rwgs  of  water  gas  shift  reaction  (mol  m-3  s-1 ) 

rp  mean  pore  radius  of  electrode  (pm) 

R  universal  gas  constant  (8.3145  J  mol-1  K-1 ) 

Si  entropy  of  species  i  (i  represents  H20,  H2,  and  02) 

Sm  source  term  in  continuity  equation  (kg  m-3  s-1 ) 

SXl  Sy  source  terms  in  momentum  equations  (kg  m-2  s-2) 
ST  source  terms  in  energy  equations  ( W  m-3 ) 

Ssp  source  terms  in  species  equations  (kgm-3  s-1 ) 

T  operating  temperature  (K) 

U  velocity  in x  direction  (ms-1) 

Uin  gas  velocity  at  the  SOFC  inlet  (m  s-1 ) 

V  SOFC  operating  potential  (V);  velocity  in y  direction 
(ms-1) 

X  molar  fraction  of  species  i 

Y  mass  fraction  of  species  i 

s  electrode  porosity 

§  electrode  tortuosity 

<7ij  mean  characteristic  length  of  species  i  and  j 

£2d  dimensionless  diffusion  collision  integral 

p  density  of  the  gas  mixture  (kg  m-3 ) 

fi  viscosity  of  gas  mixture  (kg  m-1  s-1 ) 

Baa, a  activation  overpotential  at  anode  (V) 

Bact,c  activation  overpotential  at  cathode  (V) 
rj  ohmic  Ohmic  overpotential  of  the  electrolyte  (V) 


2.0  is  usually  recommended  in  order  to  avoid  carbon  formation 
[28,42]. 

In  this  study,  the  performance  of  a  planar  SOFC  running  on 
FI20/CFl4  mixture  is  investigated  with  a  2D  numerical  model. 
It  is  found  that  running  on  FI20/CFl4  mixture  without  any  pre¬ 
reforming,  the  upstream  of  the  SOFC  works  in  electrolysis  mode 
thus  consuming  electricity  generated  in  the  downstream,  at  an 
operating  potential  of  0.8  V.  This  electrolytic  effect  in  SOFC  is  due 
to  the  non-uniform  electrolyte  Nernst  potential  and  uniform  oper¬ 
ating  potential  along  the  SOFC  channel.  This  finding  is  contrary  to 
our  common  understanding  that  electrochemical  reactions  should 
cease  under  the  open-circuit  voltage  conditions.  The  results  of  the 


Table  1 

Parameters  used  in  simulation. 


Parameter 

Value 

Operating  temperature,  T  (K) 

1173 

Operating  pressure,  P  (bar) 

1.0 

Electrode  porosity,  e 

0.4 

Electrode  tortuosity,  £ 

3.0 

Average  pore  radius,  rp  (pm) 

0.5 

Anode-supported  electrolyte 

Anode  thickness  da  (pm) 

500 

Electrolyte  thickness,  I  (pm) 

100 

Cathode  thickness,  dc  (pm) 

100 

Height  of  gas  flow  channel  (mm) 

1.0 

Length  of  the  planar  SOFC  (mm) 

20 

Thickness  of  interconnect  (mm) 

0.5 

Inlet  velocity:  Uin  (ms4) 

1.0 

Cathode  inlet  gas  molar  ratio:  02/N2 

0.21/0.79 

Anode  inlet  gas  molar  ratio:  H20/CH4 

0.667/0.333 

SOFC  operating  potential  (V) 

Thermal  conductivity  of  SOFC  component  (W  m_  1  K_1 ) 

0.8 

Anode 

11.0 

Electrolyte 

2.7 

Cathode 

6.0 

Interconnect 

1.1 

present  study  suggest  that  operating  the  SOFC  at  a  lower  potential, 
or  use  partially  pre-reformed  CH4  fuel  are  beneficial  to  eliminate 
the  electrolytic  effect  to  improve  the  electric  output  of  CH4-fueled 
SOFC. 

2.  Model  development 

In  this  study,  H20/CFl4  mixture  with  a  molar  ratio  of  2.0  is  used  in 
the  SOFC  anode  for  power  generation  via  direct  internal  reforming. 
The  computational  domain  and  working  mechanism  in  a  planar 
SOFC  is  shown  in  Fig.  1.  The  computational  domain  includes  the 
two  interconnectors,  the  fuel  gas  channel,  nickel,  yttria-stabilized 
zirconia  (Ni-YSZ)  anode,  YSZ  electrolyte,  and  YSZ-LSM  (lanthanum 
strontium  manganite)  cathode,  as  well  as  the  air  gas  channel.  Typi¬ 
cal  dimensions  of  the  gas  channel  and  the  cell  component  thickness 
are  used  in  the  modeling  study  and  summarized  in  Table  1. 

During  operation,  FI20/CFl4  mixture  and  air  are  supplied  to  the 
anode  and  cathode  channels,  respectively.  In  the  anode,  CFI4  under¬ 
goes  DIR  to  produce  FI2  and  CO  (Eq.  (1 )): 

CH4  +  H2O^CO  +  3H2  (1) 

Due  to  a  high  steam  to  carbon  ratio  (2.0  in  the  present  study), 
WGS  reaction  occurs  in  the  porous  anode  as  well  (Eq.  (2)): 

C0  +  H20^C02  +  H2  (2) 

The  H2  fuel  produced  from  DIR  and  WGS  reaction  is  transported 
through  the  porous  anode  of  the  SOFC  to  the  triple-phase-boundary 
(TPB)  at  the  electrolyte-electrode  interface,  where  FI2  molecules 
react  with  oxygen  ions  (O2-)  to  produce  H20  and  electrons  (Eq. 
(3)).  The  electrons  produced  are  transported  through  external 
circuit  to  the  cathode  to  produce  useful  power.  In  the  cathode, 
02  molecules  are  transported  through  the  porous  cathode  to  the 
electrolyte-cathode  interface  and  react  with  electrons  (coming 
from  anode)  to  produce  oxygen  ions  (Eq.  (4)),  which  are  subse¬ 
quently  transported  through  dense  electrolyte  to  the  anode  side: 

H2  +  02“  -*  H20  +  2e~  (3) 

0.5O2  +  2e~  ->■  02“  (4) 

In  addition  to  FI2,  both  CO  and  CFI4  may  be  electrochemically 
oxidized  in  the  porous  anode.  However,  since  the  electrochemical 
oxidation  rates  of  CO  and  CH4  are  very  small  compared  with  DIR 
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Fig.  1.  Computational  domain  of  a  planar  SOFC  running  on  H20/CH4  mixture. 


and  WGS  reactions,  electrochemical  oxidation  of  CO  and  CH4  can 
be  safely  neglected  [54-56].  In  addition,  chemical  reaction  of  C02 
with  CH4  is  assumed  to  be  negligible  in  the  present  study. 

Based  on  the  working  mechanism,  a  2D  numerical  model 
is  developed  to  capture  the  coupled  transport  and  chemi¬ 
cal/electrochemical  reactions  in  a  planar  SOFC.  The  model  consists 
of  an  electrochemical  model,  a  chemical  model  and  a  computa¬ 
tional  fluid  dynamics  (CFD)  model.  The  details  of  the  sub-models 
are  described  in  the  subsequent  sections. 

2.1.  Electrochemical  model 


2.2.  Chemical  model 

The  chemical  model  is  developed  to  calculate  the  rates  of  DIR 
and  WGS  reaction  and  the  resulted  heat  sink/generation.  In  the 
literature,  the  rates  of  DIR  and  WGS  reaction  can  be  determined 
approximately  as  [58-61  ]: 

Rdir  =  (1031.40238  -  1.9844T 

+(9.55556E  -  4)T2 )(1 .0E  -  8)PCh4Ph2o  (8) 

-(1.4)(1.0E-20)Pco(Ph2)3 


The  electrochemical  model  is  used  to  calculate  the  current  den¬ 
sity  (J)-v oltage  (V)  relationship.  It  is  assumed  that  the  operating 
potential  is  constant  while  the  current  density  varies  along  the 
gas  flow  channel.  This  is  valid  as  interconnector  with  high  elec¬ 
trical  conductivity  is  placed  along  the  gas  channel  to  collect  the 
current  produced  in  an  SOFC.  In  operation,  the  J-V  relationship  can 
be  described  by 

V  =  E  —  rjact,a  —  hact,c  ~  P  ohmic  (5) 

pL  (p'  )0'5’ 

H2  °2J 


where  E  is  the  equilibrium  potential  (Nernst  potential),  depend¬ 
ing  on  operating  temperature  and  gas  composition;  T  is  the 
local  temperature  (I<);  rjohmic  is  the  ohmic  overpotential  of 
the  electrolyte  (YSZ);  rjact,a  and  rjacttC  are  the  activation  over¬ 
potentials  at  the  anode  and  cathode,  respectively;  R  is  the 
universal  gas  constant  (8.3145  J  mol-1  K-1);  and  F  is  the  Fara¬ 
day  constant  (96,485  C  mol-1).  P1  is  the  partial  pressure  at  the 
electrode-electrolyte  interface,  which  means  that  the  concen¬ 
tration  overpotentials  are  included  in  the  Nernst  potential  (E). 
The  ohmic  overpotential  can  be  easily  calculated  by  the  Ohm’s 
law.  The  activation  overpotentials  are  usually  calculated  with  the 
Butler-Volmer  equations.  However,  as  the  activation  overpotential 
and  current  density  usually  follow  a  linear  relationship  according 
to  experimental  observations  [57],  the  activation  overpotentials  are 
calculated  as 

Pact  =JRact  (7) 

where  Ract  is  the  resistivity  (£2m2)  of  the  electrode  due  to  elec¬ 
trochemical  reaction.  Based  on  the  previous  study  on  activation 
overpotential,  the  overall  resistivity  of  the  anode  and  cathode  of 
0.2  C2  cm2  is  used  in  the  present  study. 


RT 

E  =  1.253  -  0.0002451 6T  +  —  In 

2  F 


Kwgs  =  (0.02377 T  -  23.84836)(1.0 E  -  7)PcoPh9o 

(9) 

-(0.03262T  -  33.46016)(1. 0E  -  7)PCo2Ph2 

where  RDir  and  RWG s  are  the  rates  (mol  nrr3  s-1 )  of  DIR  and  WGS 
reaction,  respectively;  and  p  is  the  partial  pressure  (Pa). 

The  reaction  heat  for  DIR  and  WGS  reaction  can  be  determined 
by  the  enthalpy  change  of  the  reactions  [62].  Assuming  linear 
dependence  on  temperature,  the  reaction  heat  (J  mol-1 )  for  DIR  and 
WGS  reaction  can  be  calculated  as 

Hdir  =  -(206,  205.5  +  19.51757)  (10) 

HWgs  =  45,  063  -  10.287  (11) 

The  negative  and  positive  signs  in  Eqs.  (10)  and  (11)  indicate  that 
the  DIR  and  WGS  reaction  are  endothermic  and  exothermic,  respec¬ 
tively. 


2.3.  Computational  fluid  dynamics  model 


The  fundamental  transport  phenomena  in  an  SOFC  include  fluid 
flow,  heat  transfer  and  mass  transfer.  In  an  SOFC,  laminar  flow 
conditions  are  usually  applied  due  to  a  small  Reynolds  number. 
The  transport  phenomena  in  SOFCs  are  governed  by  the  conser¬ 
vation  laws  for  mass,  momentum,  energy,  and  species,  which  are 
summarized  below  [63,64]: 


d(pu)  d(Pv) 

dx  dy 


(12) 


d(pUU)  djpvu ) 

dx  +  dy 


dp  d_(  du\  d_ 
dx  +  dx  y  dx  J  +  dy 


+  sx  (13) 


d(puv )  d(pw) 

dx  +  dy 


dp  d_ 
dy  +  dx 


(14) 


d(pcpUT)  d(pcpVT) 
dx  dy 


(15) 
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Table  2 

Parameters  used  in  calculating  the  effective  diffusion  coefficients  [65]. 


CO 

co2 

h2 

o2 

ch4 

n2 

H20 

(X  i 

3.69 

3.941 

2.827 

3.467 

3.758 

3.798 

2.641 

silk 

91.7 

195.2 

59.7 

106.7 

148.6 

71.4 

809.1 

dx 


d(pUYi)  ,  d(pVYi)  _  9  ( nn  dYj)  ,  d  fnn  dYt 
+  dy  ~  dx  1' P  l’m  9x  J  +  9y  lP,'m  9y 


+ 


»sp 


(16) 


where  U  and  V are  the  velocity  components  in  x  and y  directions;  p 
and  ii  are  the  gas  density  and  viscosity  of  the  gas  mixture,  respec¬ 
tively,  which  depends  on  local  temperature  and  gas  composition;  k 
and  cp  are  the  thermal  conductivity  and  heat  capacity,  respectively; 
Y;  denotes  the  mass  fraction  of  species  i,  which  can  be  calculated  as 


>< 

II 

(17) 

N 

m = yxM 

(18) 

In  the  porous  electrodes,  effective  heat  conductivity  and  heat 
capacity  are  used  [66]: 


k  =  ekf  +  (  1  -  s)ks 

Cp=SCpJ  +  (  1  -S)Cp,s 


(27) 

(28) 


i=l 


where  kf  and  ks  are  the  heat  conductivity  (W  m-1 1<-1 )  of  the  fluid 
and  solid,  respectively;  cpj  and  cPiS  are  the  heat  capacity  (J  kg-1 1<-1 ) 
of  the  fluid  and  solid,  respectively. 

The  source  term  5  in  continuity  equation  (Eq.  (12))  represents 
the  mass  change  due  to  electrochemical  reactions.  As  electrochemi¬ 
cal  reactions  are  assumed  to  occur  only  at  the  electrode-electrolyte 
interface,  the  source  term  is  non-zero  at  the  electrode-electrolyte 
interface  and  zero  in  other  regions.  At  the  anode-electrolyte  inter¬ 
face,  the  source  term  for  continuity  equation  can  be  written  as 
[66]: 


where  Xj  and  Mz  are  the  molar  fraction  and  molecular  weight  of 
species  i,  respectively. 

Deff  is  the  effective  diffusion  coefficient  of  species  i  in  gas 
mixture  (both  anode  and  cathode)  and  can  be  evaluated  by  Eqs. 
(19)— (21 ): 


D 


1 

w 


1 


1 


Di,m  = 


Di,m  +  Dik 


1  -Xf 


n  P 

ui,k  =  ~y 


Ej^iVDu 

2  rv  ISRT 


jrMj 


(19) 


(20) 


(21) 


where  ij/e  is  the  ratio  of  tortuosity  to  porosity  of  porous  electrodes; 
and  rp  is  the  radius  of  pores.  Dp  is  the  binary  diffusion  coefficient  of 
species  i  and  j,  which  can  be  determined  as 


Dij  = 


0.0026 T1-5 

P\/Wjalj^D 


M  -  =  _ _ _ 

lJ  (1/Mi)  +  (1  /Mj) 


(22) 


(23) 


where  cr^  is  the  mean  characteristic  length  of  species  i  andj;  £2D  is  a 
dimensionless  diffusion  collision  integral,  which  can  be  calculated 
as 


°i  +  G\ 


„  1.06036 

~  - n  icci - f 


(24) 


0.193 


■  +  ■ 


1.03587 


+  ■ 


1.76474 


0.1561  exp(0.47635r)  exp(1.52996r)  3.8941  lr 

(25) 


r  = 


1<pT 

£i  j 


(26) 


where  kb  =  1 .38066  x  1 0-23  (J  K_1 )  is  the  Boltzmann’s  constant.  The 
values  of  0[  and  £fj  can  be  used  in  the  present  study  are  summarized 
in  Table  2  [65]. 
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Fig.  2.  Electrolytic  effect  in  SOFC  running  on  H20/CH4  mixture  at  an  operating 
potential  of  0.8  V— (a)  distribution  of  current  density;  (b)  distribution  of  electrolyte 
Nernst  potential. 
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where  Aact  is  the  active  area  for  electrochemical  reaction  at  the 
anode-electrolyte  interface  and  Vc  is  the  size  of  control  volume. 
Ay  is  the  width  of  the  control  volume  in  y  direction  at  the 
anode-electrolyte  interface.  At  the  cathode-electrolyte  interface, 
this  source  term  is  calculated  as 


0.004  -i 


Sm  —  — 


JMo2 
4  F  Ay 


(30) 


where  the  negative  sign  above  means  oxygen  is  electrochemically 
consumed. 

The  Darcy’s  law  is  used  as  source  terms  in  momentum  equations 
(Eqs.  (13)  and  (14))  so  that  the  momentum  equations  are  valid  for 
both  the  gas  channels  and  the  porous  electrodes  [64]: 


S*  = 


IXU 


Sy-  B, 


(31) 


(32) 


The  source  term  ( W  nrr 3 )  in  energy  equation  ( Eq.  ( 1 1 ))  includes : 
( 1 )  heat  generation  due  to  electrochemical  entropy  change  and  irre¬ 
versible  overpotentials;  (2)  heat  energy  demand  for  DIR  (Eq.  (1)); 
and  heat  generation  due  to  WGS  reaction  (Eq.  (2)).  In  the  porous 
anode,  both  DIR  and  WGS  contributes  to  the  source  term: 


Sp  =  RdirHdir  +  ^wgs^wgs 


(33) 
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Fig.  3.  Distribution  of  chemical  reaction  rates  in  SOFC  at  an  operating  potential  of 
0.8  V— (a)  direct  internal  reforming  of  CH4;  (b)  water  gas  shift  reaction. 
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Fig.  4.  Distribution  of  gas  composition  in  SOFC— (a)  H2  molar  fraction;  (b)  H20  molar 
fraction;  (c)  CH4  molar  fraction. 


In  the  dense  electrolyte,  the  source  term  includes  the  irre¬ 
versible  loss  through  entropy  change  and  activation  losses  via 
electrochemical  reaction  (at  the  electrode-electrolyte  interface), 
as  well  as  ohmic  overpotential  of  the  electrolyte.  These  losses  are 
assumed  to  evenly  distributed  in  the  electrolyte  and  thus  the  source 
term  can  be  written  as 


JT  AS  Jr)t 

Sj=  2fT  +  T 


(34) 
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Fig.  5.  Distribution  of  temperature  in  SOFC  running  on  H20/CH4  mixture. 

where  L  is  the  thickness  of  electrolyte;  AS  is  the  entropy  change  of 
the  electrochemical  reaction  (J  K-1  mol-1 )  and  can  be  calculated  as 


AS  =  Sh2o  -  0.5So  -  Sh2 


(35) 


where  Sh2o»5o2  andSH2  are  entropy  ofH20,02,  and  H2,  respectively, 
which  can  be  found  in  Ref.  [62].  The  negative  AS  indicates  that  heat 
is  released  from  the  electrochemical  reaction. 
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Fig.  7.  Distribution  of  temperature  in  SOFC  running  on  H20/CH4  mixture  at  an 
operating  potential  of  0.4  V. 


Tjt  (V)  is  the  total  overpotential  losses,  which  can  be  calculated 

(36) 


as 

ijt  =  V  —  E 


The  source  terms  (Ssp)  in  species  equations  (Eq.  (16))  can  be 
calculated  in  the  way  similar  to  the  source  term  for  continuity  equa- 
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Fig.  6.  .  Electrochemical  performance  of  SOFC  at  an  operating  potential  of  0.4  V— (a) 
distribution  of  current  density;  (b)  distribution  of  electrolyte  Nernst  potential. 


Fig.  8.  .  Distribution  of  chemical  reaction  rates  in  SOFC  at  an  operating  potential  of 
0.4  V— (a)  direct  internal  reforming  of  CH4;  (b)  water  gas  shift  reaction. 
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tion  (Eq.  (12)).  However,  in  the  anode,  the  DIR  and  WGS  reaction 
should  be  included. 

3.  Numerical  methodology 

The  governing  equations  are  solved  with  the  finite  volume 
method  [67].  The  boundary  conditions  and  the  detailed  calcula¬ 
tion  procedures  can  be  found  in  the  previous  publication  [64].  The 
electrochemical  model  and  the  chemical  model  are  linked  with  the 
CFD  model  through  the  source  terms  in  Eqs.  (12)-(16).  In  the  itera¬ 
tion,  the  electrochemical  model  and  the  chemical  model  are  solved 
to  calculate  the  current  density  and  chemical  reaction  rates,  which 
are  used  to  determine  the  source  terms  for  the  CFD  model.  Subse¬ 
quently,  the  CFD  model  is  solved  to  update  the  temperature,  gas 
composition,  fluid  velocity  in  SOFC,  which  are  used  to  solve  the 
electrochemical  model  and  the  chemical  model.  Computation  is 
repeated  until  convergence  is  obtained.  The  in-house  CFD  code  is 
written  in  FORTRAN. 

4.  Results  and  discussion 

The  electrochemical  model,  chemical  model  and  the  CFD  model 
have  been  validated  respectively  in  the  previous  publications,  by 
comparing  the  modeling  results  with  data  from  the  literature 
[61,64].  Numerical  simulations  have  been  carried  out  to  ensure 
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Fig.  9.  Comparison  of  SOFC  electrochemical  performance  at  inlet  gas  velocity  of 
0.5  m  s-1  and  5.0  m  s-1  —(a)  distribution  of  current  density;  (b)  distribution  of  elec¬ 
trolyte  Nernst  potential. 
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Fig.  10.  Distribution  of  H2  molar  fraction  in  SOFC  running  on  H20/CH4  mixture— (a) 
inlet  gas  velocity:  0.5  m  s-1 ;  (b)  inlet  gas  velocity:  5.0  m  s-1 . 


grid  independence  of  the  results.  The  dimensions  and  typical  struc¬ 
tural/operating  parameters  used  are  summarized  in  Table  1.  As 
mentioned  in  Section  1,  a  carbon-steam  ratio  of  2.0  is  used  in  the 
present  study  in  order  to  avoid  carbon  deposition.  In  terms  of  molar 
fraction,  0.667  and  0.333  are  used  for  H20  and  CH4  at  the  anode 
inlet. 

4.1.  Electrolytic  effect  in  SOFC  running  on  H20/CH4  mixture 

One  significant  finding  of  the  present  study  is  that  the  current 
density  of  SOFC  running  on  H20/CH4  mixture  increases  sharply 
from  a  negative  value  near  the  SOFC  inlet  to  be  on  the  order  of 
1000  Am-2  in  the  downstream  of  SOFC,  at  an  operating  potential 
of  0.8  V  (Fig.  2a).  The  negative  current  density  near  the  inlet  indi¬ 
cates  that  the  upstream  of  the  SOFC  works  in  the  electrolysis  mode. 
The  positive  current  density  in  the  downstream  of  SOFC  indicates 
that  the  downstream  of  SOFC  works  in  fuel  cell  mode  and  provides 
the  electrical  energy  needed  for  H20  electrolysis  near  the  inlet,  as 
the  electrical  conductivity  of  the  current  collector  (interconnec¬ 
tor)  placed  on  the  entire  SOFC  channel  is  very  high.  As  can  be  seen 
from  Fig.  2b,  the  electrolyte  Nernst  potential  increases  from  below 
0.8  V  near  the  inlet  to  be  slightly  higher  than  0.9  at  the  outlet  of 
the  SOFC.  The  low  Nernst  potential  is  mainly  because  the  partial 
pressure  ratio  of  H2  to  H20  is  low  near  the  inlet  (Eq.  (6)).  The  non- 
uniform  local  Nernst  potential  and  the  uniform  operating  potential 
cause  the  electrolytic  effect  found  in  the  present  study,  since  the 
operating  potential  (0.8  V)  is  higher  than  the  local  Nernst  potential 
in  the  upstream  but  lower  than  the  local  Nernst  potential  in  the 
downstream  of  SOFC  running  on  H20/CH4  mixture.  This  finding  for 
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SOFC  has  not  been  reported  before  because  almost  all  modeling 
studies  on  CH4  fueled  SOFC  assumes  that  part  of  CFI4  is  exter¬ 
nally  reformed  (typically  about  30%  pre-reformed)  before  feeding 
to  the  anode  of  SOFC  (i.e.  [37-41  ]).  The  present  study  also  indicates 
that  at  open-circuit  voltage  condition,  the  electricity  produced  in 
the  downstream  of  SOFC  is  exactly  used  up  through  electrolysis  in 
the  upstream  of  SOFC.  It  is  contrary  to  conventional  understanding 
that  electrochemical  reactions  should  cease  at  open-circuit  voltage. 
However,  this  finding  is  similar  to  the  experimental  observations 
in  Prof.  Zhao’s  work  on  direct  methanol  fuel  cells  (DMFCs)  [68-71  ]. 
In  their  works,  electrolytic  hydrogen  evolution  was  observed  in  a 
DMFC  at  open-circuit  voltage  [68-71  ].  The  reason  is  that  part  of  the 
DMFC  works  in  fuel  cell  mode,  producing  electricity,  while  part  of 
DMFC  works  in  electrolytic  mode,  consuming  electricity  and  pro¬ 
ducing  H2  gas,  which  experimentally  supports  the  finding  of  the 
present  study. 

In  order  to  better  understand  the  electrolytic  effect  and  the 
working  mechanisms,  the  distributions  of  DIR  and  WGS  reac¬ 
tion  as  well  as  gas  composition  in  the  SOFC  are  studied.  In  the 
porous  anode,  the  rate  of  DIR  decreases  considerably  from  about 
388  mol  m-3  s-1  near  the  inlet  to  be  less  than  100  mol  m-3  s-1  in 
about  7  mm  downstream  from  the  inlet  (Fig.  3a).  In  the  further 
downstream,  the  rate  of  DIR  does  not  vary  much  along  the  channel. 
For  comparison,  a  locally  high  rate  (over  80  mol  nrr3  s-1)  of  WGS 
reaction  is  observed  in  the  porous  anode,  and  in  most  of  the  anode 


0.004  n 


0.003 


>0.002 


I 


H20  Molar  Fraction 


0.001 


Operating  temperature:  1173K 
Operating  potential:  0.8V 
Inlet  gas  velocity:  0.5  m.s*1 


KO/CH  molar  ratio  at  the  inlet  of  anode:  2.0 

2  4 


0.000 

0.1 


I 


0.427 

0.457 

0.487 

0.518 

0.548 

0.578 

0.609 

0.639 

0.669 


000  0.005  0.010  0.015  0.020 

Location  along  the  main  flow  stream,  x(m) 

(a) 


0.004 


0.003 

? 

>0.002 


0.001 


0.000 

0.000  0.005  0.010  0.015  0.020 

Location  along  the  main  flow  stream,  x(m) 

(b) 


0.635 


Operating  temperature:  1173K 
Operating  potential:  0.8V  0  646 

Inlet  gas  velocity:  5.0  m.s'1  mh  0.657 

HJD/CH  molar  ratio  at  the  inlet  of  anode:  2.0 
2  4  ■  0.668 


Fig.  11.  Distribution  of  H20  molar  fraction  in  SOFC  running  on  H20/CH4  mixture— (a) 
inlet  gas  velocity:  0.5  m  s-1 ;  (b)  inlet  gas  velocity:  5.0  m  s-1 . 
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Fig.  12.  Distribution  of  direct  internal  reforming  rate  SOFC  running  on  H20/CH4 
mixture— (a)  inlet  gas  velocity:  0.5  m  s-1 ;  (b)  inlet  gas  velocity:  5.0  m  s-1 . 


the  rate  of  WGS  reaction  is  about  30-50  mol  nrr3  s-1  (Fig.  3b).  The 
calculated  reaction  rates  are  on  the  same  order  but  a  little  higher 
than  the  data  by  Lehnert  et  al.  [58].  The  higher  reaction  rate  for  DIR 
and  WGS  may  be  caused  by  a  higher  temperature  and  a  higher  H20 
molar  fraction  used  in  the  present  study. 

Due  to  the  electrolytic  effect  near  the  inlet,  the  DIR  and  WGS 
reaction,  H2  molar  fraction  increases  considerably  from  0.0  at  the 
inlet  to  be  higher  than  0.23,  despite  of  electrochemical  reduction 
of  H2  in  most  part  of  the  SOFC  (Fig.  4a).  The  molar  fraction  of  CH4 
decreases  along  the  SOFC  channel  due  to  DIR,  as  shown  in  Fig.  4b. 
The  decrease  in  SOFC  temperature  along  the  channel  indicates  that 
the  heat  sink  due  to  endothermic  internal  reforming  of  CH4  is  higher 
than  heat  generation  due  to  irreversible  losses  and  the  WGS  reac¬ 
tion  (Fig.  5).  The  temperature  difference  between  the  inlet  and 
outlet  is  more  than  73 1<  in  the  present  simulation.  The  temperature 
distribution  in  the  present  study  is  different  from  an  SOFC  running 
on  partially  reformed  CH4  gas  mixture,  i.e.  work  done  by  Aguiar 
et  al.  [39].  In  an  SOFC  fueled  with  a  gas  mixture  of  CH4,  H20,  H2, 
CO,  C02  (partially  pre-reformed),  the  temperature  decreases  near 
the  inlet,  reaches  the  minimum,  and  increases  in  the  downstream 
[39].  This  means  that  if  CH4  is  partially  (i.e.  30%)  pre-reformed,  the 
heat  generation  due  to  WGS  reaction  and  electrochemical  reac¬ 
tion  can  exceed  the  heat  demand  for  DIR  in  the  downstream  of 
SOFC. 
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Fig.  13.  Distribution  of  water  gas  shift  reaction  rate  in  SOFC  running  on  H20/CH4 
mixture— (a)  inlet  gas  velocity:  0.5  m  s-1 ;  (b)  inlet  gas  velocity:  5.0  m  s-1 . 


4.2.  Effect  of  operating  potential 

Since  the  electrolytic  effect  in  an  SOFC  running  on  H20/CH4 
mixture  is  caused  by  the  higher  operating  potential  than  the  local 
Nernst  potential,  this  effect  may  be  eliminated  by  operating  the 
SOFC  at  a  lower  potential.  As  shown  in  Fig.  6a  and  b,  the  local  cur¬ 
rent  density  increases  from  a  positive  value  along  the  SOFC  channel, 
indicating  that  the  electrolytic  effect  is  removed  since  the  operat¬ 
ing  potential  (0.4  V)  is  lower  than  the  local  Nernst  potential.  As  the 
current  density  of  SOFC  at  an  operating  potential  of  0.4  V  (Fig.  2a) 
is  significantly  higher  than  at  0.8  V  (Fig.  6a),  heat  generation  from 
electrochemical  reactions  increases,  which  results  in  a  higher  tem¬ 
perature  in  the  downstream  of  SOFC  at  an  operating  potential  of 
0.4  V.  As  can  be  seen  from  Fig.  7,  the  temperature  at  the  outlet  of 
SOFC  is  increased  from  1 1 00  K  at  0.8  V  to  1 1 33  K  at  0.4  V.  As  a  result, 
the  rates  of  DIR  and  WGS  reaction  in  the  downstream  of  the  SOFC 
are  slightly  increased  when  the  operating  potential  is  decreased 
from  0.8  V  to  0.4  V  (Fig.  8a  and  b). 

4.3.  Effect  of  inlet  gas  velocity 


This  phenomenon  is  attributed  to  the  fact  that  the  gas  composition 
has  smaller  variation  at  a  higher  inlet  gas  velocity  than  at  a  smaller 
inlet  gas  velocity,  as  more  gas  is  supplied  to  the  SOFC.  As  can  be  seen 
from  Fig.  1 0a,  the  FI2  molar  fraction  increases  considerably  from  0.0 
(inlet  of  SOFC)  to  about  0.3  (outlet  of  SOFC)  at  an  inlet  gas  velocity 
of  0.5  m  s-1 .  For  comparison,  the  H2  molar  fraction  only  increases 
to  about  0.1  at  the  outlet  of  SOFC  at  an  inlet  gas  velocity  of  5.0  m  s-1 
(Fig.  10b).  Similarly,  the  FI20  molar  fraction  in  the  downstream  of 
SOFC  is  found  higher  at  an  inlet  gas  velocity  of  5.0ms_1  than  at 
an  inlet  velocity  of  0.5  ms-1  (Fig.  11).  The  higher  ratio  of  FI2/H20 
(molar  fraction)  at  a  higher  velocity  leads  to  a  higher  electrolyte 
Nernst  potential  as  well  as  a  higher  current  density,  as  shown  in 
Fig.  9.  In  addition,  the  smaller  variation  in  gas  composition  at  a 
higher  inlet  gas  velocity  results  in  higher  DIR  (Fig.  12)  and  WGS 
reaction  (Fig.  13)  in  the  downstream  than  at  a  smaller  inlet  gas 
velocity,  since  more  reactants  are  supplied  at  a  higher  velocity.  For 
example,  at  an  inlet  gas  velocity  of  5.0  m  s-1 ,  the  rates  of  DIR  and 
WGS  reaction  near  the  outlet  of  SOFC  are  about  100  mol  m-3  s-1 
and  40  mol  m-3  s-1,  respectively  (Figs.  12b  and  13b).  At  a  smaller 
inlet  gas  velocity  (0.5  m  s-1 ),  these  two  reaction  rates  near  the  SOFC 
outlet  are  about  50  mol  m-3  s-1  and  20  mol  m-3  s-1,  respectively 
(Figs.  12a  and  13a). 

5.  Conclusion 

A  numerical  model  is  developed  to  study  the  performance  of 
SOFC  running  on  FI20/CFl4  mixture  considering  DIR  and  WGS  reac¬ 
tion  in  the  porous  anode.  It  is  found  that  at  an  operating  potential  of 
0.8  V,  or  at  the  open-circuit  voltage,  the  upstream  of  SOFC  works  in 
electrolysis  mode  while  the  downstream  of  SOFC  works  in  fuel  cell 
mode.  The  electrolytic  effect  in  SOFC  is  caused  by  the  non-uniform 
local  Nernst  potential  and  uniform  operating  potential  along  the 
SOFC  channels.  At  the  open-circuit  voltage,  both  H20  electrolysis 
and  fuel  cell  reactions  take  place  in  the  SOFC  running  on  H20/CH4 
mixture  and  the  electricity  generated  from  the  fuel  cell  part  is 
exactly  used  up  by  the  electrolysis  part.  This  finding  is  contrary  to 
conventional  understanding  that  electrochemical  reactions  should 
cease  at  open-circuit  voltages. 

The  electrolytic  effect  in  the  SOFC  can  be  eliminated  by  operating 
the  SOFC  at  a  sufficiently  low  potential— lower  than  the  local  Nernst 
potential.  At  a  lower  operating  potential,  the  current  density  is  con¬ 
siderably  increased,  which  in  turn  results  in  higher  heat  generation 
through  electrochemical  reaction  and  thus  a  higher  temperature 
in  the  downstream  of  SOFC.  An  increase  in  inlet  gas  velocity  from 
0.5  m  s-1  to  5.0  m  s-1  does  not  contribute  to  the  reduction  of  elec¬ 
trolytic  effect.  Instead,  the  electrical  output  of  SOFC  is  reduced  at 
a  higher  inlet  gas  velocity,  as  both  the  average  current  density  and 
Nernst  potential  are  decreased. 

This  study  provides  better  understanding  of  the  CH4-fueled 
SOFC  with  internal  reforming  and  indicates  that  partial  pre¬ 
reforming  of  CH4  is  beneficial  for  SOFC  considering  the  electrolytic 
effect. 
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As  another  important  operating  parameter,  the  inlet  gas  velocity 
is  varied  to  examine  its  effect  on  SOFC  performance.  It  is  found 
that  an  increase  in  inlet  gas  velocity  from  0.5  m  s-1  to  5.0  m  s-1  has 
negligible  effect  on  the  current  density  and  Nernst  potential  near 
the  inlet  but  significantly  decrease  both  the  current  density  and 
Nernst  potential  in  the  downstream  of  the  SOFC  (Fig.  9a  and  b). 
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